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The sputtering process has been investigated by simulating 

the sputtering of single-crystal copper with 1-7 keV argon. A 

digital computer was used to build the crystal, bombard it, and 

move crystal atoms. Four mechanisms were observed which 

cause surface atoms to sputter. An atom is sputtered when 

(1) it is squeezed out of the surface, (2) it is scooped out 

when another atom strikes its inner hemisphere, (3) it is ejected 

when an atom passes behind it, and (i+) it is knocked out by a 

second layer atom which is moving outward. Nearly all sputtered 

atoms were surface atoms. Second and third layer atoms were 

/ 

sputtered only for ion energies greater than 5 keV. They 'were 
sputtered by mechanisms similar to the surface atom mechanisms. 
"Silsbee chains^' were observed to be directed into the crystal, 
and momentum focising was observed to cause sputtering only 
when it occurred in close packed, surface rows . Outward 
directed chains were not observed. Sputtering deposit patterns, 
sputtering ratios, and sputtered atom energy distributions were 



obtained for (100), (110), and (111) surfaces. All data 
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1 . INTRODUCTION. 



Material is lost from a substance undergoing ionic bombard- 
ment. This is sputtering . Many authors have tried to define 
sputtering using terms such as "( atoms) ... ejected or knocked 
out", "emission", "ionic erosion", "disintegration", and "break- 
down." A reader who blithely accepts one of these definitions 
in the course of his reading may conclude: This definition 

describes the mechanism of sputtering. This is an illogical, if 
not erroneous, conclusion; the definition of sputtering must be 
derived from the mechanism, not the mechanism from the 
definition. The sputtering process has been investigated 
utilizing one of two general mechanisms, ejection or emission . 

Ejection implies an immediate or direct reaction release of an 

» 

atom from a crystal; emission, a delayed or indirect reaction 
process which results in the release of an atom. 

In 1923, Kingdon and Langmuir bombarded thoriated 
tungsten with various ions in a glow discharge tube. This was 
a special case of sputtering since the thin surface film of 
thorium on a tungsten substrate was sputtered rather than the 
tungsten itself. The results of this experiment, that the 
sputtering ratio, atoms removed per incident ion, increased 
with increased ion mass and increased ion energy, qualitatively 
suggested an ejection mechanism. A few years later, Von 

( P ) 

Hippel and Blechschmidt' ' proposed a theory which described 
sputtering as an evaporation of surface atoms, an emission 
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mechanism. Earlier, Von Hlppel^^^ had found by spectroscopic 
means that at least some sputtered atoms were in an excited 
state. The sputtering theory showed that atoms in the region 
of impact could rapidly acquire thermal energy if the kinetic 
energy of incident ions was converted to thermal energy at the 
target surface. If an atom acquired a sufficient amount of 
thermal energy, it would then evaporate from the surface, some 
atoms evaporating while in an excited state. This theory was 
improved by Townes^^^ in 1944> who calculated a sputtering 
flux and the number of atoms evaporated per incident ion. 
Keywell^^^ in 1954 » used neutron diffusion theory to approxi- 
mate atomic interactions within the crystal, a new approach to 
the theoretical investigation of sputtering. Direct application 
of statistical methods to sputtering was made by Harrison^ 
who envisioned the Interaction of two distribution functions 
(the crystal lattice and the Ion beam). These models, based on 
statistical methods, implicitly accept ejection type mechanisms, 

* P 

One of the most Important contributions to the study of 
sputtering was made by Wehner^'^^ in 1953* In the first 
sentence of his paper, Wehner stated, "The most widely 
accepted sputtering theory is the evaporation theory...", but 
the deposit patterns of single crystal sputtering which he 
obtained showed pronounced, high density areas or "rpota". 

There was now strong evidence for a momentum transfer 
process, and further development of the evaporation theory 



14 



ceased. Shortly after Wehner’s findings were reported, 
Henschke^®) proposed a theory of sputtering based solely on 

t 

classical collision theory, treating normal and oblique incidence 
sputtering separately. Oblique incidence sputtering could be 
explained by two body collision processes. The incident ion 
penetrated the surface layer with little or no interaction, was 
relected outward by second layer atoms, and ejected a surface 
atom by striking it on its inside hemisphere. His concept of 
normal incidence sputtering required many-body collisions in 
which the ion was eventually reflected outward to sputter 
surface atoms. This theory was plausible for oblique incidence 

i. 

sputtering, however, the case for normal incidence sputtering 
required that the ion be reflected inside the crystal. (For an 
ion more massive than the target this requirement cannot be 
met . ) 

Silsbee noted that in any discussion of momentum trans- 

fer effects in sputtering, the geometry of the crystal should 
be considered; successive collisions in a crystalline structure 
might be influenced by the structure itself. His calculations 
showed: in a two-body collision, an atom has a departure angle 
with respect to an axis specified by a close packed direction. 
In subsequent collisions, the departure angle , will decrease 
if the energy is low enough, and the collision sequence occurs 
along a close packed row of atoms. This concept of "momentum 
focusing", "Silsbee chains", or "focusons" was considered by 
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many investigators, excepting Wehner^^^), to be a fully satis- 
factory explanation of the spots in deposit patterns. But the 
deposit pattern is only one of many observable features of 
single crystal sputtering. 

The sputtering ratio is another important characteristic 
of the sputtering process. If a theory or model is to explain 
sputtering, both patterns and sputtering ratios must be 
explained qualitatively and quantitatively. Alm^n and Bruce 
(AB) measured sputtering ratios of a variety of metals 

using N, Ne, Ar , Kr , and Xe ions over a 5—65 keV range of 
bombardment energies. (Although polycrystalline specimens 
■were sputtered, the qualitative results are the same as those 
obtained for single crystal sputtering.) They noted that the 
sputtering ratio was an increasing function of - bombardment 
energy provided the ion mass was greater than the target mass. 
Sputtering ratios showing a maximum, or of a slowly varying 
nature, were observed in all cases in which the ion was lighter. 
Evidence such as this coincides with any gross conception of an 
ejection mechanism. 

Single crystal copper was sputtered by argon at inter- 
mediate energies by Magnuson and Carlson and South- 

ern, Willis, and Robinson ( SWR ) ' . MC measured sputtering 

ratios (using 1-10 keV argon) for the (ill), (100) and (110) 
surfaces and found that sputtering ratios decreased in this 
order of surfaces. SWR sputtered single crystal copper with 
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1-5 keV argon and measured sputtering ratios, but of greater 
significance was the quality of their deposit patterns which 
clearly showed the presence of only certain spots: 

(ill) surface 3-(H0)> l-(lll) but not seen due to 

beam aperture. 

(110) surface l-(llO), 2-(l00) 

(IGO) surface 4-(H0), I-(IOO) 

The work of AB , MC, and SWR indicated that the sputter- 
ing ratio is at least a function of bombardment energy (or 
momentum) and of mass ratios. However, consideration of 
spot patterns indicates that the sputtering process itself is 
also highly dependent on the gross and/or surface geometry of 
the crystal and on ion penetration depth. 

Ion penetration of solids has been investigated experiment- 
ally and by computer simulation. Piercy, McCargo , Brown, 
and Davies' ^ investigating channeling of various heavy ions 
in monocrystalline aluminum, found that the ion penetration 
distance increased in the order (111), (100), and (llO) beam 
orientations. This is the same order as the sputtering ratio 
decrease found by MC. A number of studies of ion penetra- 
tion have also been made using computer simulations (Oen, 
Holmes, and Robinson^ Robinson and Oen^ , Harrison, 

Leeds, and Gay( ^ . 

The most apparent correlation between the results of ion 
penetration and sputtering studies is the concept of trans - 
parency , defined by Fluit, Rok and Kistemaker ^ ^ . Both 
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experimental and simulation investigations of ion penetration 
indicate that channels exist in certain crystallographic directions 
of a crystal. Results of normal incidence sputtering experi- 
ments indicates that sputtering ratios are lower for surfaces 
in which these channels are parallel to the beam. One might 
initially suspect that sputtering ratios will be low if ions and 
target atoms are confined to these channels; momentum is 
directed into the crystal rather than laterally, reversal of 
momentum to produce sputtering not occuring. Accordingly, 
some incident ions should see a transparent surface rather 
than a wall of atoms. 

The obvious complexities of the sputtering process may be 
investigated separately by digital computer simulation. This 
of course involves the selection of a gross mechanism ( ejection 
or emission) and a model (hardsphere or otherwise) . If 
factors such as the interatomic potential function and its 
parameters are known, then the use of the right model in the 
computer may be expected to simulate the actual sputtering 
process. Explanation of deposit pattern characteristics and 
sputtering ratios will necessarily follow from an exact 
simulation . 



18 



2 . 



PURPOSE OF THE INVESTIGATION. 



The belief that the formation of. spots in sputtering deposit 

; 

patterns is solely a consequence of momentum focusing was 
widely accepted for a number of years. The importance of 
momentum focusing, not only in spot formation, but in the 
entire sputtering process is now believed to be much less than 
originally thought This investigation was undertaken to 

explore the single crystal sputtering process using a digital 
computer to simulate a copper crystal which is being bombarded 
with a beam of argon ions. Primary effort was directed 
towards; 

a. Finding the mechanisms which result in atoms being 
deposited in the spots or high density regions of the 
deposit pattern. 

b. Determining the correlation between crystal structure 
and sputtering ratios and patterns. 

c. Investigating the surface binding energy and its effect 
on deposit patterns. 
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3 . MODEL , 



A. A space lattice is established "whose sites represent 

the equilibrium positions of copper atoms in a face centered 

cubic crystal. Atomic spacing is that determined by x-ray 

o 

crys"callographic studies (for copper, a = 3*615 A). The 
potential function, with parameters for copper-copper inter- 
action, is the Born-Mayer type Gibson Number Two . The 
interaction between argon and copper is similarly described 
using a Born-Mayer potential, V(r) = e , but the numerical 

values of the parameters A and for this function are those 
determined by Harrison, Carlston, and Magnuson (HCM)(2^) 
from a study of secondary electron emission. The range of both 
potential functions is eroded at rQ, one half the nearest neigh- 
bor distance such that the potential and force go to zero for 
atoms whose centers are farther apart than 2r^ . Use of this 
eroded form of the potential functions permits establishment 
of the crystal in a stable, static state without use of an 
attractive potential. This equilibrium state is destroyed only 
when an atom moves from its lattice site. 

The lattice as a whole is restricted by the current computer 
program to perfect form. Defects such as interstitial atoms 
cannot be used, however, irregularities in the surface layer are 
introduced by using two types of variation from a perfect 
surface; vacancy and stub conditions. The perfect or regular 
surface is identical to a full lattice plane. In the vacancy 
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configuration the atom which would normally be hit first by 
the ion is removed from the perfect surface. The stub 
surface has an atom added on the perfect surface in a stable 
position and adjacent to the target atom. Additional atoms 
can also be added on the surface or removed from the perfect 
surface to provide a variety of random irregularities. 

B. Each ion of the beam is approximated by a single, 
neutral argon atom whose velocity vector intersects the plane 
of surface atoms at an impact point . (The term ion is used 
throughout this thesis for the incident particle to avoid con- 
fusion with target atoms . ) 

Successive runs are made using different impact 
points to simulate use of the entire beam area. The use of 
this simple model for the ion beam is subject to four 
conditions ; 

1. The results obtained using one random set of impact 
points are essentially the same as those obtained 
using any other random set of impact points. 

2. The majority of beam ions are neutralized before 
reaching some arbitrary boundary which defines the 
surface. 

3 . Prior to neutralization, the path of a beam ion is 
not appreciably affected by the surface potential. 

If. A crystal region is in an equilibrium state each time 
a neutralized ion impacts in that region. 
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Ea ch of these conditions has been satisfied either 
in the model itself or by calculation. A full discussion is 
contained in Appendix A. 
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4. DYNAMICS. 



The force functions for copper-copper and argon-copper 
interactions are derived from the respective potential functions. 
Atom position and velocity are then determined using Newton's 
Second Law of Motion. Normally, one would determine these 
values by numerical integration methods over small time inter- 
vals, however, this process uses too much computer time. 
Special methods approximating integration have been developed 
which are used to minimize computer run time yet maintain a 
good approximation to the integration. 

A. Gay developed an iterative method of solution which 
is similar to that used by Gibson, Goland, Milgram, and 
Vineyard^ Both Gay's and Gibson's methods use the 

technique of replacing differential quantities with finite 
differences to approximate the equation of motion, F/m = dv/dt 
t','/ Vt by F/m = Av/ £^t. Using this equation and the relation- 
ships Csx = v/it and v = (v + v^)/2, equations - (3-1) and (3-2) 
are obtained. 

Av = V - Vq = j^F (x^)/mj ^t (3-1) 

Ax = X - x^ = |^[f (x^)/m^ At/2 + (3-2) 

The known values F(x^), v^, x^^, and ^^t in these equations 
are usually associated with a time t = t^; x and v> the 
unknown values, are to be evaluated using these equations at 
time t = T(t^, At)'. 

Gibson's technique associates the value of x^ with a time 
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t “ t , but V is associat.ed with a time t = t - ^t/2. 
ox 0 o ov o 

Current values of x and v are alternately computed at intervals 
At using equations (3-1) ar d (3-2). 



The advantage of this technique over that of evaluating both 
variables at time t “ t^ + ^ t is that a continuous smoothing 
of the values of x and v occurs with each successive computa- 
tion. 

Gay recognized the need for some sort of smoothing 
procedure but considered that an averaging process for equation 
parameters rather chan computed results would result in a 



better approximation. His method replaces the evaluated 
force function F(xq) in these two equations with an arbitrary 
force function f ( x) which is linear in some interval to be 
determined. Consider the following two step cycle: 

Move an atom from x^^ to x^ using computed values of 
F(xq) and v^ to solve equations (3-1) and (3-2). Compute 
F(xj) and average this value with F(x^). Now, move the 
atom from X] to xg using the averaged force, F(xq, x.j^) 
and Vq to solve the equations using the same constant value 
of At. L-et f now be defined in the interval (x^, •^^)> such 
that f(xo) >-< F(Xq, ) . Equations (3-1) and (3-2) may 
now be written: 



v(t^ + £^t/2) ■'^'o “ ^^2) + ^F(xQ)/m"^ t 

x(tQ +At) - x(t^) + 




V ( t -f At) V ( t ) 
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x(t +^t) = x(t) ( XQ)/rnj^t/2 + v(t)j ^ t 

The values of v and x are determinable in any interval in 
which f approximates ^ to some arbitrary degree, therefore, 
in every complete cycle or timestep a new interval and new 
function f are defined. 

A determination of which of these two methods is better 
must be based on some standard or condition. In many cases 
the accuracy of numerical solutions such as these may be 
compared to the exact value one obtains by integration. But, 
there are no truly exact results in a simulation since a physical 
model is in the computer rather than the computer being used 
as a means of solving equations. However, if the model is 
physical, physical laws must be satisfied, and the energy gain 

or loss due to the mathematical approximations must be small. 

( 31 ) 

This condition of energy balance has been used by Harrison'"^ ' 
who found that, in simulations involving atomic interactions such 
as channeling or sputtering, use of Gay's method resulted in a 
better energy balance. 

B. The value of these approximation tecliniques may be 
substantially reduced by injudicious selection of a numerical 
value for At. T oo large a value invalidates the approximation, 
but too small a value increases computer running time. The 
program used in this simulation of sputtering incorporates not 
only Gay's approximation method but also a procedure which he 
developed for automatic adjustment of At. Its value is 
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calculated using the relationship Tt = Sr/v where ^t is replaced 
here by ^t to indicate that its value varies with each timestep. 
The velocity, v is the magnitude of the velocity of the most, 
energetic atom in that timestep, and Sr is used here as an 
input parameter, the timestep multiplier . This parameter is 
assigned the computer variable name DTI and is defined: the 
maximum displacement of any atom in any timestep shall be the 
numerical value of DTI in appropriate units of length. Gay 
considered the value of ^t, consequently that of DTI, to be a 
function of ion energy, impact point, and start point, but 
he was unable to find the relationship. Trial and error methods 
were used to find optimum values, and reasonably good results 
were obtained using these values. Johnson^^^) reported that 
the energy balance was maintained to within 3% for all his runs. 

The difficulties in choosing a value for DTI have been 
largely eliminated in this study. Some problem areas such as 
start point and impact point dependence were avoided by always 
positioning the ion tangent to the first atom It would hit. 

This, at worst, will result in a constant error for which 
adjustments may be made (this has not been necessary). The 
positioning procedure is described in Appendix B . The depend- 
ence on ion energy is Inherent to the DTI computation process 
and is easily explained once thS clutter of Impact point and 
start point dependence has been cleared away. Use of Vhe 
DTI process excludes all variable dependence except the Inverse 
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proportionality to velocity, but the interval in which the force 
function is considered linear is exactly the same width as the 
numerical value of DTI. Accordingly, a straightforward 
analysis is made to show the interdependence. 

The force function, F(r), is expanded about an arbitrary 
point r* using a Taylor’s series for a function of one variable: 
F(r)=F(r'l<) +1)F 

r r'l= 2!'^„2 k!”3r^r* 



It 

( r-r* ) 



Since F(r) is derived from the potential function V(r) = gA+Br, 

A ..T - k 

o '' ^ , F(r) is necessarily of exponential form and » F = B F 



r^ 



S' , and the series is rewritten in the following form: 



F(r) = 



|^F(r>:=) jl + B(r-r*)|| + ^F(r>l<) B^ (r-r*)^1^ (3-3) 

The quantity (r-r^^) is now defined as a variable displacement 
^r , and equation (3-3) may be written: F(r) = f( S r) + g( 5 r) 
where f ( S r) and g( S r) are the terms in braces in the force 
function expansion. The function f (cTr) is assumed to be the 
linear function used in Gay's approximation method; the function 
f ( r) is to approximate F(r). Thus, it is required that the 
ratio'ig( ^r)/f{S r) be some fraction less than one. 

Let an assumption be made that the ratio is much less 
than one, but that the actual value is the fractional deviation ’ 
from linearity, a variable whose value is to be specified. Then, 
since B is a negative constant (from the potential function), 
the truncation error for the series is easily found, and the 
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quantity B d r is known for every specified value of fractional 
deviation. But B is constant, and once the fractional devia- 
tion is specified, r is fixed. Consequently, specified to the 
fractional deviation from linearity, there is a function f(r) 
which is a valid approximation to the function F(r) in every 
interval (r, r+ir) where ^r is fixed. The quantity ^ r is 
thus the value to be used for DTI, and t = ^r/v ~ DTl/v. 

It is now assumed that Lt will always be the maximum 
possible value which minimizes computer running time without 
introducing errors due to the non-linearity of the true force 
function. One would expect therefore, that the energy 
balance would be maintained to the same order as the linearity 
of the force function. Further consideration will show that 
the energy balance is maintained to a considerably better 
degree than that anticipated. First, since DTI is the maximum 
displacement of the most energetic atom , essentially all other 
atoms will be displaced a distance less than DTI. Thus, the 
interval in which the force function f(r) is used is smaller, 
and f is an even better approximation to F_. Second, for a 
given fractional deviation from linearity, the value of ^ r is 
total displacement, not x> X.> or* £ component displacement 5 
and the velocity used in the determination £^t - DTI/v for 
each timestep is the total velocity v, not the component 
velocity vj . The result of these considerations may be 
summarized by a calculation for the maximum displacement of 
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an atom along the i component of its displacement vector: 

^x. = v-C^t, but V. 4 V ^ ^'^^/'''’max* follows 

that itiX; = Vj (DTI,/vj^^^) 4 DTI. 

Tnese analyses appear to provide a tidy solution to the 
time optimization and energy balance problem. This would be 
true except for the complications introduced by the eroded 
form of the potential function. Consider a situation in which 
one atom is separated from another by a distance ( 2r^ + d) 
where d is some distance less than the total distance D that 
the atom will move in the next times tep. When the atom is 
moved the distance Dj as shown in figure 1, the force on the 
atom at its new position is the same force it would have if 
it had moved only (D-d) units. Since the force is zero for 
r>ErQ, the model assumes that the atom had a velocity v = 
(D-d)/ t (which is smaller than its actual velocity v = D/^ t , 
and a smaller kinetic energy is computed using this smaller 
velocity. This difficulty was originally foreseen by Gay who 
included corrections to the force calculations for such situa- 
tions. Both Levy ^30) a.nd Johnson^ improved the methods 

of corrections, and recent improvements in the force calcula- 
tions have been made by Harrison^^^. 

The continuous improvement of the dynamics section of 
the program has resulted in a model in which not only has the 
computer running time been minimized, but of greatest signifi- 
cance, for nearly all of the 700—800 combinations of' surfaces. 
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ion energies, and impact points used in this study, the energy 
balance was maintained to within 1 % . 
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5. RESULTS . 



Reports of sputtering studies often include ad hoc formula- 
tions of mechanisms which are used to explain deposit patterns, 
sputtering ratios, and sputtered atom energy distributions. 

The sputtering process is discussed in this thesis in terms of 
mechanisms which have been observed in the simulation to sputter 
atoms . The observance of these mechanisms is an advantage 
peculiar to a simulation. Each crystal atom (and the ion) must 
be identified by number, at least for purposes of computer 
calculations. Since every atom is identified, its complete track 
can be plotted and labeled. The tracks of selected atoms can 
then be superimposed to show a complete set of interactions. 

The mechanisms which were observed to cause an atom to 
sputter are considered the prime observable quantity of the 
simulation. 

The results of the simulation are discussed for each crystal 
surface. The determination that an atom having energy _E is 
sputtered, was made using a probability-of-sputtering function, 
P(E) = 1 - e-E/Eb v>Jhere E^ is the assumed binding energy. 

The value of P is compared with a random number R having a 
value between 0 and 1. If P ^ R , the atom is considered 
sputtered. (This method of selection has been compared with 
one using a-^step fimction in which all atoms with E ^ E^ are 
assumed sputtered. A selection method using another probabil- 
ity-of -sputtering function, P(E) = 1 - e ' •“ b ' has also been 
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tested. No differences were observed in the sputterlai 



patterns based on the atom leleetlon method.) 

Mechanisms observed to cause sputtering are discussed 
first for each surface. Deposit patterns, sputtf^ring ratios, 
and sputtered atom energy distributions are then discussed. 
The deposit pattern^ shown are superpositions of those 
obtained by sputtering the crystal separately with regular, 
vacancy, and stub surface conditions. This practice follows 
that used by Johnson^ who found that pattern features 
were difficult to recognize when viewing Separate patterns 
from each surface condition. The energy distributions are 
also superpositions of the three surface conditions. The use 
of superposition Is considered to most realistically simulate 
the condition of the crystal surface at various times during 

f 

sputtering. There are strong arguments however, for using 
only the regular surface condition. Sputtering ratios are 
not averaged (which would correspond to the superposition 
of other data) since a probability factor for the existence 
of each surface condition would need to be determined. In- 
stead, the simulation sputtering ratios for vacancy and stub 
surface conditions are Included and discussed briefly only for 

t 

completeness, (The (110) surface was sputtered only at i 

e 

and 3 k©V using all three surface conditions.) 

Crystallographic nomenclature Is used generlcaliy In the 
discussion of spots, lieference to an (hkl) spot does not 



imply ^hkl^ ejection . The ^ hkl^ directions are with respect to 
a right handed coordinate system in which the crystal is 
described. The (hkl) surface is contained in the x-z plane; 
the 21 direction is into the drystal. The three crystals used 
in the simulation are shown in figures 2-4; all contain I 50 
atoms. This crystal size has been found sufficient to contain 
nearly all energetic collisions for ion energies up to 7 keV. 

A. The (111) surface. 

The (111) surface was sputtered normally with argon at 
1-5, 7, 10, 20, and 40 keV. (Although the potential function 
is considered valid only for ion energies less than about 7 keV, 
the higher energy runs were used to search for additional 
mechanisms.) Four mechanisms were found. These were 
especially evident at the lower (1-5 keV ) ion energies. Three 
of these are classified surface mechanisms since only surface 
atoms were found to participate in the sputtering event. The 
fourth mechanisms is a deep mechanism in which an inward 
moving atom is reflected and then initiates a sputtering event. 
In all cases, the dominant mechanisms were found to be 
surface mechanisms. 

Except for head-on collisions, the ion is scattered by the 
target (atom 2) with a component of momentum parallel to 
the surface as well as a much larger pei’pendicular component 
into the crystal. The effect of this combination of momenta 
is to drive a nearest neighbor (n. n. ) , for example atom 6 in 
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figure 5, into an apparent channel formed by atoms 7, 

8, and 9* The channel is apparent since it terminates 
abruptly thus causing atom 6 to drive atoms 7 and 8 into the 
crystal and atom 9 outward. Since atom 6 passes nearly 
directly behind atom 9» the impulse is more normal to the 
surface than parallel to it. This is termed a mole mechanism 
since one atom burrows between two layers parallel to the 
surface to sputter an atom in the outermost layer. Atoms 
sputtered by this mechanism are almost always a next nearest 
neighbor (n.n.n.) to the target and located in or near the 
sextant defined by the impact area. They usually sputter 
with greater than 10 ev. The target always receives the 
majority of energy transferred by the ion, but it is normally 
driven into the (111) trigonal array of atoms directly behind 
it. These three atoms act as buffers and dissipate the 
target's energy and momentum into the crystal . The target 
is reflected, but it does not retain sufficient energy to 
sputter a surface atom nor sputter itself. 

The second and third surface mechanisms occur when the 
impact parameter is about one third of an atom radius or 
less. In both mechanisms a n.n., atom 6, is struck by the 
ion such that the parallel and perpendiciilar momentum compon- 
ent magnitudes which it acquires are of the same oi-der. If 
the perpendicular component is greater, atom 8 scoops a n.n.n., 
atom 9, relatively high up in the surface plaiie of atoms and 
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into their outer hemispheres. These atoms are densely packed 



in the surface and act as a rigid reflector for atom 9* It is 
found that atoms sputtered by the scoop mechanism have 
ejection angles nearer the normal than would be found in the 
absence of reflection in the surface. If the perpendicular 
momentum component of atom 6 is nearly the same or less 
than the parallel component, it squeezes atom 9 against its 
surface neighbors. The squeezing causes the surface plane to 
warp and atom 9> 10, or 11 is sputtered. Atoms sputtered 
by either the scoop or squeeze mechanism usually sputter with 
less than 10 ev. Atoms sputtered by combinations of the 
three surface mechanisms are found to sputter with higher 
energies than when sputtered by a single mechanism. 

The fourth mechanism requires that an atom be reflected. 
The atom which is most often apt to be reflected is the 
target. If the impact parameter is greater than one half 
the atom radius, the target is driven towards an edge of the 
trigonal array behind it rather than into it. It will penetrate 
this edge if it has sufficient energy, but it will lose most of the 
energy during the penetration. The reflection occurs from 
third layer atoms and the target atom’s reversed momentum is 
transferred to a second layer atom. The second layer atom 
will then sputter one of two surface atoms at energies upwards 
of 5 ev. At lower ion energies the target atom cannot 
penetrate the array edge and is reflected from the second 



35 



layer; it may or may not sputter a surface atom, but it 
rarely sputters itself . 

When the crystal was sputtered at 7 keV, the ion was 
found to penetrate well into the second layer. This is a 
significant increase in ion penetration depth since 5 keV ions 
are found to penetrate only to the order of half an atom radius. 
The effect of deeper ion penetration is to cause a lateral 
compression of the second layer with accompanying warping. 
Second layer atoms which are squeezed outward will sputter a 
surface atom but are rarely found to sputter themselves. At 
the higher bombardment energies tested (10, 20, 1+0 keV ) , the 
number of second layer atoms which sputtered increased nearly 
proportionally to the increase in ion energy. The sputtering 
process at these higher energies appears to be predominantly 
by the three surface mechanisms described for lower energy 
sputtering. Mechanisms peculiar to high energy sputtering have 
not been observed although this may be due to the limited 
energy range for model validity. 

All (111) simulation deposit patters showed the character- 
istic features of three (110) spots and a central (111) spot. 

The 2 keV pattern definitely showed the presence of three 
additional spots which were found in the regions in which 
streaking occurred at other energies. These patterns are 
shown in figures 6-21 . The appearance of a hexagonal pattern 
at 2 keV was a surprising result, but it is not a unique 
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occurrence in sputtering studies. Anderson and Wehner^^^^ 
found a hexagonal pattern for (111) copper sputtered by- 
mercury. The expected trigonal pattern was found for 
energies up to 4 OO ev, a hexagonal pattern appeared at 4 OO 
ev, and it disappeared as the ion energy increased. Very 
recently, Robinson and Southern^ have found additional 
spots near ^ 114 ^ positions for ( 111 ) gold sputtered with 
4 keV argon. 

The pattern features vary slowly with binding energy 
provided its value is restricted to the range 2. 50-3 ♦ 50 ev. ' 
The upper limit is the sublimation energy for copper (33)> 
above it the patterns appear to deteriorate. This was most 
evident at a 3 keV bombardment energy. A set of patterns 
at this bombardment energy is shown in figures 8 -I 4 for a 
binding energy range of 1. 50— 5*50 ev. The pattern deteriora 
tion is most evident in the region of the (ill) central spot. 
Figures 11 and 12 show the transition from a well defined 
central spot at 3*00 ev to complete deterioration of the 
central spot at 4*00 ev. 

The half -intensity width of the (110) spots is estimated 
from pattern data and from numerical data to be about 11 
degrees. The width of the central spot is an unreliable 
datum and was not determined. 

Sputtered atoms were always found to be surface layer 
atoms for sputtering at 5 keV or less. The atoms which 
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•were sputtered most frequently appeared to naturally group by 
energies into two categories; atoms with perpendicular energy 
greater than 3 but less than 10 ev and atoms with perpendicular 
energy greater than 10 ev. The higher energy atoms were 
found to be equally distributed between the three (110) spots 
and the central spot. They were also found in the narrow 
sectors defining the streaks and at distances from pattern 
center corresponding to the distances of the (110) spots. 

This distance was found to be 0.7 units, the distance measured 
for the spots in the patterns of SWR . The lower energy 
atoms were usually found in the central spot which explained 
its sensitivity to binding energy. 

The correlation of sputtering frequency and crystal 
location of sputtered atoms was found using f requency-location 
diagrams. These are shown in figures 22-29 regular 

surface sputtering. Similar diagrams have been made for 
sputtering the vacancy and stub surfaces. It was found that 
a 1 keV ion does not cause frequent sputtering of n.n.n'. 's or 
n.n.'s in the sextant defined by the impact area. However, 
at 2 and 3 keV, n.n.n. ’s are found to be sputtered most 
frequently, and at higher energies, both n.n.n. 's and n.n.’s 

are frequently sputtered. The three-fold relationship between 

( 

the frequency of sputtering, location of the atom with respect 
to the target, and the ion energy may be clarified by consider- 
ing the sputtering mechanisms. 
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At 1 keV and an impact parameter of about three fourths 
an atom radius, the ion is found to penetrate less than half 
an atom radius. At energies up to 5 keV the penetration is 
only slightly greater. It was found that, although the target 
is always driven into the crystal, the n.n. which initiates the 
majority of sputtering may not be driven far enough into the 
crystal to provide the scoop mechanism. If it is not driven 
inward a sufficient distance to scoop the n.n.n., its energy 
will be propagated in the surface along close packed rows 
which originate at the n.n.n. This is not a ^^focuson.^ ^ At an 
early point in the propagation the surface will have warped 
sufficiently to cause an atom to sputter. The sputtering is 

a result of a squeeze mechanism but not one directly involving 

i 

the n.n. As the ion energy is increased the scoop and mole 
mechanisms are more apt to occur although the squeeze 
mechanism is found to occur at all ion energies. The transi- 
tions from dominance of one mechanism to another with 
increasing ion energies may be inferred from the profiles shown 
in figure 31 for atoms which sputter most frequently. 

A relatively constant sputtering frequency is observed for 
atoms 102 and 132 , but the percentage of high and low energy 
sputs varies for atom 102 which is a n.n.n. The percentage 
does not vary for atom 132 which is neither a n.n.n. nor a 
n.n. Atom 102 sputters most frequently at low ion energies 
by the squeeze mechanisms, at higher ion energies by the scoop 
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and mole mechanisms. These relationships have been confirmed 
by the atom track displays. They also have shown that atom 
132 is sputtered by atom 102 through a squeeze mechanism 
alone at low ion energies but in conjunction with the mole 
mechanism at higher ion energies. The profile shown for atom 
86 is the result of simple reflection from second layer atoms. 
Atoms 87 and 101, which are not in the impact sextanr but 
are sputtered frequently, are sputtered by the squeeze 
mechanism . 

Figure 32 shows profiles of atoms sputtered by the deep 
mechanism. At low ion energies the target reflects at an 
oblique angle from second layer atoms rather than penetrating 
the layer, and it enters one of the apparent ^112^ channels. 
Atoms 25 and/or 55 are then sputtered directly by the target. 
(Although this is similar to a mole mechanism it is a distinct 
mechanism since the target must be reflected in this case but 
not in the former case.) At higher ion energies the target 
will penetrate the second layer and is reflected from third 
layer atoms. It does not channel, but it causes a second 
layer atom to sputter surface atoms. Accordingly, atom 55 is 
sputtered by the deep mechanism more often than atom 25 is 
sputtered by a channeled target. Atom 40 or 70 is 
sputtered by a squeeze mechanism, usually in conjunction with 
the sputtering of atoms 25 and 55* 

These results for the sputtering of a regular surface 
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have been found to be generally applicable to sputtering vacancy 
and stub surfaces. Equal numbers of atoms from the three 
surfaces are found in the spots and in the streak regions. 

It is also found that these atoms are usually the same from 
all three surface conditions. They differ only in energy since 
it is the ion penetration distance which indirectly determines 
the magnitude and direction of momentum transferred to 
surface atoms. The vacancy surface sees an ion which 
penetrates deeper; this is analogous to an ion of greater energy. 
In the stub condition, the ion does not penetrate the surface; 
it transfers energy to the stub atom. The stub then assumes 
the role of the ion impinging on the surface, but the stub 
appears as an ion with less energy. 

Figure 33 shows the sputtering ratios for each surface 
condition as a function of ion energy; assumed binding energy 
is a parameter. A statistical variation of +20^ of the 
sputtering ratio has been assumed. This would normally be an 
unreasonably large deviation for laboratory results, but it is 
considered conservative for the numerical results from a 
simulation in which many parameters are unknown. The actual 
value of sputtering ratios obtained for the regular surface is 
generally low by 2-3 atoms/ion for a binding energy of 3-5 ev. 
Further investigation showed that, in a few cases, surface 
atoms 3 planes below the impact point would be sputtered with 
about i+-5 ev. The average increase in the (ill) sputtering 
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ratio over the 1-7 keV range is estimated to be no greater 
than 2.0 atoms/ion for a binding energy range of 1.5 to 3-5 
ev. This range is shown by the dotted lines in figure 32 a. 
Vacancy and stub sputtering ratios show erratic behavior over 
the energy range studied. 

Energy distributions of sputtered atoms are shown in 
figures 35-39 • The preponderance of atoms with energies less 
than 5 ev is caused by the lack of an intrinsic surface binding 
energy. The sputtering selection process considers all atoms 
having _E greater than 1 ev. Thus the P(0.50) energy, using 
an assumed binding energy of 3*5 ev, is 2.42 ev. This is 
especially evident for the 1 and 5 keV distributions. If this 
region is disregarded, the peaks which occur at 5-8 ev for 
all ion energies are assumed to be the maximum for each 
distribution. Secondary maxima appear at 16 ev for 1 keV 
sputtering, 16 ev for 2 keV, and at 14 ev for 3 keV sputter- 
ing. At higher ion energies, a secondary maximum may be 
present at 47 ev for 5 keV sputtering. The number of atoms 
sputtered in the simulation is too small to make a definite 
statement regarding maxima. Differences are usually measured 
by one, at most, two atoms. 

B. The (100) Surface. 

The (100) surface was sputtered normally with 1, 3» 5» 
and 7 keV argon. Mechanisms observed to cause sputtering 
were found to be nearly identical to the surface mechanisms 
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discussed for (111) sputtering. A. mechanism directly 



comparable to the deep mechanism was not observed. Instead, 
it was found that the scoop and mole mechanisms are enhanced 
by the presence of the (llO) channels parallel to the surface. 
Surface atoms which are driven into these channels do not 
have to burrow between the first and second layer in order 
to sputter nearby surface atoms. Additionally, it was found 
that these atoms will frequently sputter themselves by 
reflection from second layer atoms. The squeeze mechanism 
was observed to be the most effective sputtering mechanism 
at 3 and 5 keV. At 1 and 7 keV ion energies, variations of the 
scoop and mole mechanisms were dominant. The low-high and 
middle ion energy dependence observed for the dominant 
mechanism is directly related to the shape of the impulse by 
which the tion transfers energy to a n.n. 

At low ion energies, the impulse is sufficiently broad that 
the n.n. is directed into the edge of the square of atoms 
behind it. The n.n. will enter the (110) channel without 
immediately sputtering a surface atom. When the impulse is 
narrow, the peak force is not necessarily greater. The target 
has received most of the energy given up by the ion, and the 
ion-n.n. impact parameter is larger. The n.n. will then scoop 
or squeeze a surface atom causing it to sputter. At high ion 
energies (7 keV ) the ion was found to penetrate at least into 
the second layer where it was deflected towards but not into 
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a (110) channel. This was observed to cause second layer atoms 



to be scooped up and sputter surface atoms. Second layer atoms 
are also squeezed by their neighbors and may sputter through 
the vacancy left by a sputtered surface atom. 

Sputtering at 3 and 5 keV was found to be caused by the 
scattering of surface atoms along the surface as the ion was 
being reflected from second layer atoms. The squeeze 
mechanism was observed to sequentially sputter atoms along 
a close packed surface row. The occurrence of a mole 
mechanism was conspicuously rare; even at 5 keV the energy 
propagation was clearly restricted to propagation parallel to 

4 

the surface layer of atoms and to propagation Into the crystal. 
The only definite momentum reversals observed were for the 
target atom and the ion. 

The characteristic features of (100) sputtering deposit 
patterns are four (110) spots and a central (100) spot. The 
patterns are usually outlined by a hypocyclold-shaped haze; 
the (110) spots form the cusps of the hypocyclold. Simula- 
tion patterns were found to show these characteristic 
features quite well for 3 and 5 keV sputtering. The 3 keV 
pattern (figure 1 ^. 1 ) is a very good likeness to one at c'.b keV 
reported by SWR . The (110) spot distances from pattern 
center were found to be 1.0 TC unit which corresponds to 
the distances In the 2.5 keV experimental pattern. It was 
further determined from numerical data that the (TIO) spots 



•were actually elliptical rather than circular. This latter 
characteristic is undoubtedly the result of distortion inherent 
in flat plate collection of sputtered atoms ^ The 3 keV 

pattern was found to remain essentially unchanged as the 
binding energy was varied from 2.50 to 3*50 ev. The 5 keV 
pattern (figure i+2) was almost identical to the one at 3 keV 
only when a binding energy near 2.00 ev was used. As the 
binding energy was increased, the (110) spots became less 
well defined (figure 2+3) and spots corresponding to ^210^ 
became the most prominent feature. The intensification of 
^210^ regions was accompanied by a loss in definition of the 
hypocycloid outline. Similar results were obtained when sputter- 
ing at 7 keV. The hypocycloid outline in the pattern could 
only be observed for a binding energy near 1.50 ev. 

The deposit pattern for 1 keV sputtering (figure 40) 
showed very little similarity to those obtained from higher 
energy sputtering. Four spots corresponding to ^211^ were 
observed; four (110) spots may also be defined but they are 
extremely diffuse. (The (110) spots are best seen in the 
point plot in figure 40a.) The finding of (211) spots at 
1 keV, and the appearance of (210) spots for 5 and 7 keV 
ion energies suggested that the pattern might be rotated 45 
degrees for low ion energy sputtering. The possibility that a 
rotation occurs has been investigated but one has never 

been experimentally observed. 
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The different (100) pattern features which appeared with 
varying ion energies and binding energies were a marked contrast 
to the relative constancy observed for (111) patterns. (The 
use of a binding energy of 3*00 resulted in comparable (111) 
patterns for all ion energies.) The lack of constancy for the 
pattern simulations may be explained as a result of the sub- 
traction of a binding energy from the perpendicular component 
of a sputtered atom's energy. (This is done to simulate the 
energy lost in overcoming the surface potential.) An atom 
sputtered with small perpendicular energy, less than 15 ev for 
example, may suffer an apparently small change in its perpendic- 
ular velocity, but the direction of its velocity may be consider- 
ably altered. The average perpendicular energy of atoms found 
in spot regions is about 15-20 ev from sputtering the (100) 
face. It is about 20-25 ev for (111) sputtering for which 
variations of up to 1 ev binding energy have had little effect on 
pattern features. This explanation does not, however, explain 
why the pattern at 3 keV is valid for a range of binding 
energies. A consideration of the atoms sputtered may clarify 
but not completely explain this situation. The dissimilarities 
seen in the (100) patterns may be directly inferred from the 
frequency-location diagrams for (100) sputtering (figures 46 - 14 . 9 ) • 
When the energy propagation is restricted to the surface, 
as it is for 1 keV sputtering, the (111) surface can dissipate 
energy through four atoms in the momentum forward semi- 
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circle. The (100) has only three atoms comparably located 



to dissipate the energy. Track displays have shown that for 
1 keV sputtering, atoms 26, ll6, and 132 are sputtered by a 
combination of mole and squeeze mechanisms. The mole is 
predominant in sputtering the atoms nearest the target, the 
squeeze mechanism becoming predominant for surface atoms 
which are located further from the target. Atom 26 sputters 
less frequently than atom 116 since it is on the opposite side o'' 
of the crystal with respect to the impact area. Both of 
these atoms are found in (211) spots, but atom 116 is also 
found in (110) regions of the pattern. Atom 132 is found in 
(211) spots and the (100) spot; it is sputtered only by the 
squeeze mechanisms. When the bombardment energy is 3 keV, 
the deposit pattern has its expected characteristics. Correla- 
tion of high atom density regions with frequency-location data 
showed that the (110) spots are the result of sputtering of 
n.n.n.'s. Concurrently, the hypocycloid outline was found to 
be formed by: (1) the sputtering of the same atoms which 

formed (211) spots at 1 keV, and (2) the sputtering of the 
n.n.n.'s in the surface. These atoms were not densely 
deposited in the (211) regions but formed diffuse ellipses whose 
semi-minor axes were along quadrant bisectors. The appearance 
of well defined (110) spots at 3 keV is strongly suggested by 
the 3 keV frequency-location diagram, figure 47 • Atoms 70 
and 72 were observed to be sputtered with about the same 
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frequency as atoms Sb* Since atoma VO and 72 

are along a (100) axis, one might axpeot that they would be 
found in (110) apota, and this haa been obaerved. Atoms 35 
and 56 are in posltione relative to the target aueh that they 
would sputter within adjacent 45 degree sectors in quadrant-a 
II and III. These atoms are found In the hypooyelold outline. 

The absence of the hypocyelold at 1 keV and Ita presence 
at 3 , 5 and 7 keV Is readily apparent when one eonaldei'S the 
sputtering profiles In figure 50. At 1 keV Ion energy, n.n.'s 
are never sputtered, however, they are sputtered with 
Increasing frequency as the ion energy Increases. Atom 86, 
the n.n. In tlie quadrant containing the Impact area, exhibits 
this behavior of Increasing sputtering frequency. Atom 72 

A 

shows a sputtering pi'oflle which Indicates the (110) spot 
formation at ion enei'gies of 3 keV or greater. The constant 
nature of the profile for atom 102 at 3 keV and higher ton 
energies Is Indicative of the hypocyelold outline rather than 
(^'dUiy spot formation. Atoms 26, 116 and 132 show a generally 

■ e 

dscreaslrig frequency of sputtering with high Ion energy. This 
is consistent with i:he loss of (211) spots at higher Ion energies. 

'the track patterns observed for 5 aad 7 keV (100) 
sputtering showed that atoms ejected by a eguee«?e niechanletn 
are often ejected In directions opposite to those one would 
exjject solely on tlie basis of the location of the atom with 
respect to tlie Impact point. These- .oeeurrenttes wete caused 



by two factors; (l) the atom was squeezed against its 
neighbor and reflected outward and with its parallel component 
of momentum reversed rather than being strictly ejected out- 
ward. (2) Sputtering by near simultaneous squeeze and mole 
mechanisms results in ejection of an atom where the ejection 
direction is dependent only on the impulse delivered by the 
channeling atom in the mole mechanism. 

Sputtering ratios for the ( lOO )} surf acecare sho-Wmin f igtire 51 • 
The correspondence between simulation values and experimental 
values is quite good for the regular surface for a binding 
energy of 3*50 ev. The ratios for vacancy and stub surfaces 
show a closer correlation to the regular surface for (100) 
sputtering then for (111) sputtering. 

The energy distributions of sputtered atoms are shown 
in figures 52 - 55 * No specific maxima other than that near 
5-7 ev appeared to be present. 

C. The (110) Surface. 

This was the last of the three face centered cubic crystal 
surfaces to be sputtered in this simulation. It is purposely 
the last to be discussed. First, experimental deposit patterns 
from (110) sputtering show only a large central oval area for 
bombardment energies greater than a few hundred ev (10). 
Second, a unique feature of this surface is the (110) channel; 
no other (hkl) surface in a face centered cubic crystal shows 
^hkl^ channels. The sputtering mechanisms occurring for 1, 3> 
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5 and 7 keV ion energies were most readily observed for 
sputtering of this surface, and they are interrelated with the 
presence of (110) channels. The mechanisms are, again, 
identical in concept to those previously discussed for (111) 
and (100) sputtering. Before discussing these mechanisms 
separately in terms of individual crystal atoms, the features 
of the frequency-location diagrams (figures 5^-59) are 
summarized . 

Sputtering of the (110) surface differs significantly from 
the sputtering of (100) and (111) surfaces: (1) the target 
atom was the atom sputtered most frequently at all ion 
energies. (This is in agreement with results obtained by 
Levy (2) the number of sputtered atoms with high 

energies were generally those in the close packed row contain- 
ing the target atom, but with the target atom located at 
the center of the row rather than at the origin of the row. 
These findings are wholly consistent with expected results 
when one considers the sputtering mechanisms with respect to 
this particular surface. 

The frequent sputtering of the target is made possible by 
the nature of the squeeze mechanism but the mechanism does 
not itself cause the sputtering. When the target is driven 
into the crystal it squeezes the n.n. (in its row) aiid ejects 
it. Once the n.n. position is vacant, the target will travel 
through a large, potential-free area before st:riking a second 



layer atom. The target strikes this second layer atom and is 
free to reflect outward without finding a surface atom directly 
in its path. The target then transfers most of its parallel 
momentum to a n.n.n. in the row. The parallel impulse is 
then propagated down the surface row. Warping of this one 
row causes additional atoms to sputter. This sputtering 
sequence is also observed as a result of the ion striking the 
n.n. which is in the quadrant defined by the impact area. A 
rather surprising continuation of this sputtering mechanism 
was observed to occur in atom rows located both above and 
below the horizontal row containing the target. Second layer 
atoms which receive energy from the target or ion are driven 
into (111) apparent channels. They pass behind another second 
layer atom causing it to be ejected outward and strike two 
adjacent surface atoms. Sputtering is, again, initiated in a 
surface row. Sputtering profiles for the target (atom 72) 
and atoms 42, 71, 73 and 102 are shown in figure 60. Atoms 
i+2 and 102 are the n.n.'s to the target; atoms 71 and 73 are 
those which initiate sputtering in the horizontal rows containing 
them. Second layer atoms are sputtered more frequently as 
the ion energy is increased. 

Simulation deposit patterns from (110) sputtering (figures 
6 I- 64 ) tend to show more of an oval outline than a uniformly 
dense central oval region. The central region can be made 
more dense by including atoms sputtered with less than 1 ev 



51 



perpendicular energy in the probability-of-sputtering selection 



process. This did not seem justifiable for (ilO) pattern 
production since these small energy atoms were excluded from 
(100) and (111) patterns. The 1 keV point plot (figure 6la) 
showed indications of (100) spots but these cannot be clearly 
seen in the deposit pattern (figure 6lb) . Their presence was 
substantiated by numerical data. The average energy deposit 
data showed that the average energy of atoms in ^lOO'^ regions 
was nearly double that of atoms found in the central oval. 

This is a constant characteristic of non-central spot regions. 
At higher ion energies there was no direct indication of the 
presence of these spots either by numerical data or energy 
deposit data. A few atoms having high energy wei’e found in 
these regions but the area density was too small to form a 
spot in the pattern. The set of atoms found to form the 
central oval for sputtering at all ion energies was not a well 
defined group such as has been found for (lOO) and (111) 
sputtering. The atoms forming the (100) spots from 1 keV 
sputtering are those which are either n.n.'s to the target 
or n.n.’s to the atom which initiates sputtering in each 
horizontal row . 

The sputtering ratio foi* (110) sputtei’ing is shown in 
figure 65 . Agreement with experimental data from the 
numerical and curve shape aspect was poorer fo^‘ this sut'fact 
than the other two surfaces. This is apparently caused by 



3 ^ 



(110) sensitivity of the sputtering ratio to the A parameter 
in the argon-copper potential function. Levy ^^9^ obtained 
better curve shapes using A = 11.435 whereas A = 12.56 
has been used in this study. This sensitivity has not been 
observed for (100) and (111) sputtering. 

Energy distributions are shown in figures 66— 69» The 
notable feature of (110) sputtering is that the energy of 
sputtered atoms is generally higher for (110) sputtering than 
(100) or (111) sputtering. 

D. Results of Ancillary Studies. 

During the course of the simulation, it became apparent 
that certain extensions of this study should be made. It was 
particularly desirable to probe five areas: 

(1) Sputtering mechanisms peculiar to high bombard- 
ment energies ( up to 40 keV ) . 

(2) Sputtering of the (110) surface by high energy 
ions which are directed towards (110) channels. 

(3) Sputtering mechanisms for an ion heavier than 
the target. 

(4) Sputtering of the (0001) basal surface of a 
hexagonal close packed crystal. 

(5) Sputtering the (100) surface at ev ion energies. 

Investigation of the first two of these five areas were 

limited somewhat by the restrictions for the copper-copper 
and argon-copper potential function. They are usually assumed 
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valid up to about 10 keV, however, if small impact parameters 
are assumed to rarely occur, valid results may be expected. 

Areas (3) and ik) involve the use of unknown parameters. 
Accordingly, the validity of results from investigating these 
two areas cannot be assured. 

The search for mechanisms peculiar to high energy sputtering 
was made using the (ill) surface. This surface was considered 
the one most likely to show additional mechanisms since it is 
the most densely packed surface. The failure to find mechanisms 
(part A of this section) peculiar to high energy sputtering 
was not too surprising. 

The sputtering of the (110) surface by ions directed into 

(110) channels was investigated by determining the sputtering 
ratio for ions impacting at the two impact points located 
nearest the channels. The sputtering ratio curve for channel 
shots is shown in figure 70. The discontinuity between the 
low and high ion energy curves is assumed to be the result of 
breakdown of the potential functions. The important result 
is that sputtering always occurred for channel shots, even 
when 40 keV ions were used. 

The third area was examined using xenon to sputter the 

(111) surface. Argon-copper potential parameters were used 
since those for xenon-copper were not available. Runs were 
made at 3 and 5 keV . At both energies, the ion penetrated 

at least to the fourth layer of atoms; no indication of momentum 
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reversal of the ion was observed. The ion was found to be 



channeled, just before penetrating the third layer, into (110) 
channels. It initiated cascades which propagated into the 
crystal . 

The fourth area was investigated to see if the effect of 
surface geometry on the sputtering deposit patterns could be 
determined by a comparison of (111) sputtering of a face 
centered cubic crystal with sputtering the (0001) surface of 
a hexagonal close-packed crystal. Hasiguti, Hanada, and 
Yamaguchi^ ^ have sputtered zinc with 8 keV argon; the 
deposit pattern showed an outlined, equilateral hexagon with 
a central haze. An attempt was made to reproduce this 
pattern with the simulation model modified for zine-zinc 
interactions. The potential parameter was adjusted for a 
Born— Mayer type potential V(r) = Ae"*^'^^ where A = 52 

Q / 1 

( ^ keV. This relationship was determined by Ander- 

sen and Sigmund^ ^ . The argon-zinc potential parameters 
(assuming a Born-Mayer type potential) were approximated by 
using those for argon-copper interactions. Neither the 
sputtering ratio nor deposit pattern matched the reported 
results. This was not unexpected since the zinc atoms in 
the crystal must be represented as ellipsiods (with the minor 
axis in the basal plane) rather than represented as spheres. 

The (100) surface was sputtered at 100 ev to investigate 
the possibility of a 45 degree pattern rotation for low ion 



55 



energy sputtering. The number of atoms which were sputtered 
was extremely small. All sputtered atoms had energies less 
than 10 ev. It was found that the atoms sputtered most 
frequently from the regular surface were n.n.n.'s (figure 7^ ) • 
Bombarding vacancy and stub surface conditions also resulted 
in a predominance of n.n.n.’s being sputtered rather than 
sputtering of the a.toms forming the square about the target 
atom . 

E. Results Common to All Surfaces. 

Momentum focusing was observed to occur only into the 
crystal for the three copper surfaces studied. This was 
particularly evident from atom track displays of (ill) surface 
sputtering. If an atom underwent momentum reversal, it was 
always by reflection from atoms located no deeper in the 
crystal than the third layer. Even for 7 keV ion energies, 
crystal atoms located deep in the crystal were always driven 
inward. A surprisingly large numbei' of atoms with lai-ge 
energies were found to move between second and third layei's 
with their motion nearly pai-allel to the surface and with a 
small, inward-directed momentum component.. The energy 
which appeared to be delivered in impulses to atoms still at 
lattice sites, dissipated through the crystal. These atoms 
rarely acquired more than idO ev through the energy impulse 
process. Their motion was restricted since they were 
surrounded by other atoms. 
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The target atom in the (lOO) and the (110) surfaces, 
■when struck at near zero impact parameter by the ion, 
always transmitted the majority of its energy billiard-ball 
fashion in close packed rows perpendicular to the surface. 
(Atoms to which this energy was transferred escaped through 
the back face of the crystal.) There were a few cases in 
which one of these atoms would initiate another chain. The 
chain was never more than 2-3 atoms in length before the 
energy was either wholly dissipated or the chain stopped by 
divergence of the momentum to form numerous small cascades 
directed into the crystal interior. 

The unique arrangement of atoms in the first layer was 
observed to be a dominant factor in determining ejection 
directions for all surfaces. Atoms which were ejected other 
than nearly normal to the surface were always influenced by 
their neighbors. Atoms which were ejected at angles near 
45 degrees to the normal would reflect from their neighbors, 
often ending up as normally ejected atoms. This effect was 

I I 

most pronounced for (111) sputtering since the six atoms 
acted as a lens, but it was also seen for (110) sputtering 
in which a lens is formed by second layer atoms. The (100) 
surface showed a strong lens effect although one would not 
necessarily assume that it occurs. Sputtered atoms, 
originally adjacent in the (100) surface, were often observed 
to be ejected almost simultaneously. In these situations they 
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■were ejected nearly normal to the surface. This may be the 

( o cr ) 

genesis of dimers recently observed by Woodyard . 

The spot regions of che deposit; patterns from (100) and 
(111) sputtering always contained more high energy (> 10 ev) 
than low energy atoms. But the energy distribution was more 
uniform than increasing with the ion energy. It was not 
peaked at any one or group of energies. Atoms with low ion 
energies were predominantly found in the central region of the 
pattern . 
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6 . Conclusions. 



The observable quantities of sputtering appear to be 
interdependent on only the sputtering mechanisms. There is 
no evidence that these quantities are interrelated . The 
qualitative and quantitative data from the simulation indicate 
that the deposit pattern, sputtering ratiot, and energy 
distribution of sputtered atoms cannot be correlated with 
each other; one cannot predict a sputtering ratio from an 
energy distribution. But, each of these quantities can be 
cross-correlated between surfaces. 

The main features of deposit patterns appear to be 
determined only by the surface geometry of the crystal. The 
formation of ( 111 ) pattern spots is attributed primarily to 
assisted focusing by the hexagonal lens. The predominance of 
trigonal rather than hexagonal symmetry is considered a 
natural result of the brief ( 112 ) channeling observed in the 
mole mechanism. The appearance of streaks between the 
spot pairs and the appearance of a hexagonal spot pattern at 
2 keV suggests that use of a hemispherical collector in the 
simulation will show the three (llif) spots which were observed 

( QO ) 

by Southern and Robinson' . However, no distinction can 
be made between the ( 110 ) and (II4) spots on the basis of 
pattern location. If the mole mechanism is predominant at 
certain energies the ( 110 ) spots should be more intense than 
the (114) spots. If the squeeze mechanism is predominant, 
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the (110) and (lli+) spots should have more equal intensities. 

The nearly exact simulation of a (100) patrern at 3 keV 
is considered one of the best arguments chat sputrering is 
mostly a surface phenomenon. Neither the potential form 
nor its parameters for copper-copper and argon-copper inter- 
actions are known with certainty. If sputtering was a deep 
phenomenon, one might accept an argument that spots in the 
pattern could be produced without exact knowledge of the 
potential. However, the hypocycloid outline, definitely present 
in the simulation pattern, is considered the feature of the 
pattern which would be most sensitive to small variations in 
the potential if sputtering involved more than the first few 
layers. This conclusion is consistent with the apparent 
rotation of the 1 keV pattern. Atoms found in ^211^ regions 
were sputtered when the energy transferred by the ion to the 

target was as small as 1+0 ev. This compares favorably with 

o 

a calculated transfer value of 75 ev for a 0.5 A impact 
parameter (21)^ sputtering threshold enei'gy of 50 ev 

has been reported ( 32 )^ 

The inability to produce consistent (100) patterns at a 
given binding energy for 5 and 7 keV sputtering is undoubtedly 
due to the unknown condition of the surface. It would be 
unreasonable to assume that the binding energy is a decreasing 
function of ion energy only for the (100) surface. The 
disruption of the surface by the first group of incident ions 
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is a factor which cannot be ignored. The sputtering ratios 
determined for the (lOO) surface as well as for the (110) 
and (ill) surfaces, were within reasonable limits when the 
assumed binding energy is in the range 1.50 to 3.50 ev. No 
definite value of binding energy can be determined from the 
patterns and sputtering ratios unless a weighting factor is 
used for the surface condition at the time of ion impact. 

The assignment of weighting factors would be, at best, a 
guess . 

Three surprising features were observed in simulation 
deposit patterns. The first, the apparent rotation of the 
(100) pattern at 1 keV, will be further discussed. The 
rotation of the (100) pattern for low energy sputtering has 
not been observed. The results of the simulation indicate 
that such a rotation is possible. This conclusion is based on 
empirical rather than theoretical considerations. The hypocycloid 
outline found at 3 keV in the simulation appears to be formed 
by the same atoms which form (211) spots at 1 keV. 
Furthermore, no explanations have been previously proposed as 
to why the outline is a hypocycloid in experimental patterns. 

One would expect that a circular or perhaps a non-distinct 
outline would be observed when the (100) surface is sputtered. 
Accordingly, it is proposed that the hypocycloid outline is a 
result of the inability of the (100) surface to completely 
focus atoms sputtered by a squeeze mechanism for bombardment 
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energies at which the energy propagation is not contined in the 
surf ace . 

The second feature was the unique form of the bonds which 
connected the spots in the simulation pattern for (111) 
copper sputtered at 5 keV . The appearance of the similar 
outline, a spherical triangle, was observed in a pattern of (111) 
copper sputtered with 1.5 keV krypton by Yurasova and 
Bukhanov^ . The third feature was the hexagon spot pattern 
at 2 keV, a pattern seen by Anderson and Wehner^^*^^ for 
(111) copper sputtered with I4.OO ev mercury. The similarities 
observed in the simulation patterns with these anomalous 
features of the experimental patterns suggested that common 
factor other than the (111) copper surface might be present. 
The momentum ratios of i+OO ev mercury-2 keV argon and 
1. 5 keV krypton-5 keV argon are 1.0 and 0.8 respectively. 

This is not considered to definitely establish a sputtering 
correlation based only on ion momentum. The momentum ration- 
for 1.5 keV krypton and 4 and 3 keV argon are 0.9 and 1.0 
respectively. Ideally, the spherical triangle would bo observed 
in the 3 keV (ill) simulation pattern. Additionally, when 
(111) copper is sputtered at high temperatures, Uie ti'iangulai 
outline has been observed to become more pronounced . 

The absence ol' such a momentum scaling el feet couid Do 
easily shown by experimentally sputtering (ill) copper with 
ion momentum comparable to that ol 40O ev meimury . 



The general consistency of the numerical values of 



sputtering ratios, and energy distributions of sputtered atoms 
and likeness of the patterns with experimental data is a 
strong argument for the validity of a computer simulation of 
sputtering. It is remarkable that the consistency is as 
good as it is. The model uses only a repulsive potential, the 
crystal size used in the simulation is an infinitesimal portion 
of the smallest laboratory specimens, and the potential form 
and parameters are comparatively crude. It is concluded 
that the results obtained substantiate the concept of 
transparency and the occurrence of momentum focusing only 
within the surface layer. 
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APPENDIX A 



The Beam Model 

A. Impact areas and Impact points. 

Each crystal surface contains an intrinsic, plane geometric 
shape; hexagon for (111) surface, square for (100) surface, 
and rectangle for (110) surface. A volume element of the 
(hkl) surface is defined by the area of this intrinsic shape 
and a depth of some number of (hkl) planes. This volume is 
chosen so that by translation along axes of a Cartesian 
coordinate system the entire crystal may be generated. The 
smallest intrinsic area which can be chosen for each surface 
and still satisfy the translation requirement is shown in 
figure 72 . A finite number of points is symmetrically distributed 
within each of these areas to represent the infinite set of 
possible points of impact for an incident ion. These smallest 
areas are further divided into representative impact areas 
which are seen from figure 73 to be degenerate under appro- 
priate rotations and/or inversions of the coordinate axes. 

Since the sectors are degenerate, the set of impact points 
in each sector is also degenerate; only the points contained 
in one impact area need be used to represent bombardment of 
the entire area. The coordinate axes rotation and inversion 
schemes are discussed in Appendix C (Deposit Pattern 
Production ) . 

The independence of the impact point set used, with 
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respect to number and kinetic energy of atoms sputtered, was 
tested using two sets of points. No dependence was found. 
Results for the (111) surface using a 3 keV ion are shown in 
figure 74 as an example. Set 1 points are those shown in 
figure 73 > set 2 (not shown) is a set of eleven points which 
are located between the points of set 1. 

B. Neutralization of Bearn ions. 

The assumption that argon ions are neutralized prior to 

impact on copper may be inferred from results of a theoretical 

( P7 ) 

study of secondary electron emission by Harrison et. al. ' . 

Consideration of atom-atom rather than ion-atom interactions 
gave results in reasonable agreement with experimental data. 

A supporting argument for neutralization is based on 

/ O U \ 

Hagstrurn's theoi>y of Auger ejection of electrons'-^ . The 
probability of an ion being neutralized in ^ at x is; 

Pt(x,v)«a exp ^ - exp -a( x-Xj.^^ )"j -a(x-x^)'|^ 

whex’e x^^ “ (l/a) Sin (A/av) is the value of _x whore is 






maximum. 'I'he paratnetei-s A and a occur in the transition 

rate function, and jv is the velocit.y of an ion start:mg at x oo 

Hagstrurn used tungsten as an example and obtained t:ho va)iu> 

of A by empirical means, the value of a from published datai. 

Rather than atterniJt an exact proof of neutralization for argn 

on copper, Hagstrurn's result:s fox' ai’gon on tungsten lu^ve been 

used to give an order of magnitude, at worst, approxunation . 

o 

Accordingly, 1.6341x^,^2.18 A fox' 1 to 10 ke'7 ^u’gon ions on 
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copper; this is a reasonable distance from the crystal surface 
and one may assume that neutralization occurs. 

C. Ion Deflection by Surface Potential. 

As a first approximation, it is assumed that a singly 
positive-charged argon ion sees the crystal as a single, fixed 
scattering center having 2-1 positive charge. The well known 
equation for a central force-induced hyperbolic orbit is used 
in conjunction with figure 75 to determine the deflection. 
r ( 0 ) = ( L^/ mk) / ^ -1 + ( 1 + 2EL^/ mK^) cos( 

L = ion angular momentum m = ion mass 
E = ion energy K = 

= angle of closest approach 
The angular momentum L is determined at 

o 

to be L = mvb where mv‘^/2 = E and b is the impact parameter. 
The equation for r( 6 ) is put in a more convenient form by 

P 2 P 

making^ the substitution = (mvb) = 2mEb to give 

r (©) = b (2Eb/k)/|^ -1 + ^ + (2Eb/k)^^^/^ cos(Q-(9o ) 
The angle of closest approach is determined as a function 

of the product Eb by requiring the denominator to vanish for 
Once (as a function of Eb) is known, the 
deflection ratio r/b is calculated at the crystal surface. 

The deflection calculated is for an unneutralized ion, thus 
greater than that for an ion which is neutralized at some 
distance in front of the surface. Plots of r/b at the surface 
and (0 as functions of the ion kinetic energy-impact parameter 
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product are conrained in figure ?6. The numerical values are 
for an Ar — Cu system. 

A more sophisticated approximation, scattering from a 

fixed dipole of same charge sign, requires a messy integration. 

Rather than follow this line, one may expect from the nature 

of the problem that the path of the approaching ion -will undergo 

some oscillatory motion or perhaps corkscrew motion as the ion 

is influenced by the surface potential. In either case one 

would expect that the net acceleration of the ion parallel to 

the crystal surface would be no greater than that due to a 

single fixed scattering center. With these considerations, 

figure 76 is used to determine the percentage of ions which 

will be appreciably deflected. As an example of an appreciable 

ion deflection, consider r/b “ 'd.0\ from figure ?6, Eb “ 3*6 x 
o 

10“’*^ keV-A, An energy range of 1-10 keV for argon ions, 

that used in this study, corresponds to a range of impact 

o 

parameters of 0.035 - 0.0035 A. The fx’action oi surface 
area, and therefore fraction of Ions which will be appreciably 
deflected is (0 .033)^/ ( 1.^6)^ ?.? x 10~^^ for 1 keV ions i-ina 

7.7 X lO”* for 10 keV ions ( 1 . kib A is the copper atom 
effective radius in the crystal). This negligibly small iract.iv.m 
of ions cannot Influence the macro aspects of the sputtering 
model. 

D. Equilibrium £btate During Impact. 

The use of a single atom '.ipproximaLlon to the be.itu wiili 
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the crystal model described in the main text requires that the 

crystal region be in an equilibrium state at the time of an 

ion impact. It is not necessary that this equilibrium state 

be identical to previous equilibrium states since the variety 

of surface configurations available in the model provides for 

random surface conditions. It is necessary that the time 

required for the crystal region to return to an equilibrium 

state be small with respect to the arrival-time intervals of 

the ions. Satisfaction of this condition is determined by 

comparing an experimental beam flux over the area of the 

crystal face of the model used to the time required for 

completion of all energetic collisions in the model. 

Beam intensities of the order of 10(^amps/cm'^ were used 

by Magnuson and Carlston^ ^ . If a beam of this intensity 

o 

is incident on a crystal surface area of less than lO^A^ such 

as in the model, the ion flux over this area is less than 100 

ions/sec, an ion arrival-time interval of 10” seconds; all 

energetic collisions in the model are completed within about 
—12 

10 seconds. The relaxation time of the region is thus much 
smaller than the ion arrival-time intervals and the condition, 
that the crystal region be in an equilibrium state at the time 
of impact, is satisfied. 
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APPENDIX B 



Positioning the Ion 



The ion is positioned tangent to the first target atom it 
■will strike. Since this is a dynamic rather than static process 
it is not necessary that this be a stable position on the crystal 
surface. Figure 77 shows the ion at its arbitrary initial 
position and calculated final position, both with respect to the 
impact point and target atom. The initial position is a small 
distance in front of the surface, beyond the eroded potential 
range of the crystal atoms. This position is described by a 
vector r^ originating at the impact point and having direction 
parallel but opposite to the ions velocity vector. The vector 
rg, from the impact point to the target atom, is known 
since the target atom's coordinates are known; the vector r^ 
is known, and it is desired that vector r^^ have a magnitude 
equal to the distance between centers of two tangent atoms, 
2r^. The law of cosines is used to compute the magnipxde of 
r^ which lies along r^ . Accordingly, the following calculations 
are made: 

cos OC - "r . "r" 



The positive squai'o root solution is chosen to give the tarigent 
on Lhi outside hemisphere of the tai’get atom. 



1 



Kll 1^2| 
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The ion positioning is accomplished automatically for each 
run by subroutine START, which is contained in the computer 
program for the sputtering simulation. 
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APPENDIX C 



Production and Analysis of Sputtering Deposit Patterns 
A. Production. 

A sputtering deposit pattern represents the intersection 
points of atoms' velocity vectors with the surface of a 
collector plate. The sputtered atoms from the simulation 
are collected on a flat plate by determining these points of 
intersection. Each atom which exits through the crystal 
surface is initially assumed to have been sputtered and a data 
card has been prepared for each one. Data of particular 
interest are each velocity component magnitude and the kinetic 
energy perpendicular to the surface; data of secondary 
interest are' atom number, impact point used and ion kinetic 
energy. (The use of these last data will be discussed shortly.) 

It is recognized -that each sputtered atom has lost some 
energy to overcome the surface binding energy. An assumed 
value of binding energy is subtracted from the perpendicular 
kinetic energy and a new perpendicular velocity component is 
calculated. Parallel velocity components are normalized to the 
new perpendicular component to give a two dimensional cooi‘din- 
ate point. This point is the intersection of the atom’s 
velocity vector with an imaginary collector plate placed at unit 
distance from the target surface. Each point is tlien rotated 
and/or mirrored about the coordinate axes to give the inter- 
section of a velocity vector which would have resulted for aii 



Ih 



ion impacting at the corresponding impact point in each of 
the other impact areas. The impact areas are shown in 
figure 73 and the coordinate point rotation and mirroring values 
are listed in program DATASORT which is used to generate 
the points. Each point is plotted using program DATAPLOT 
with a CDCI 6 OA computer and a CalComp plotter to give a 
point plot . The dimensions of the plots are in target- to- 
collector or T-Q units since these points have been normalized 
to unit targe t-to-collector distance. The scale which has 
been used permits plotting deposit points of atoms which have 
been sputtered within an escape cone of about 63 degrees 
(57 degree cone shown in figures). This has been found 
satisfactory to contain all pattern features of interest. 

The point plots show only point patterns and therefore 
do not accurately simulate experimental patterns which are 
area density patterns. The conversion of a point plot to a 
smooth area density pattern is made by photographing the point 
plot with the camera defocused such that no single point is 
distinguishable but high and low density areas are prominent. 
Developing and printing is controlled to bring out the high 
density areas while maintaining the haze background. Loss of 
intensity in some spot regions cannot be avoided such as seen 
in figures 6b and 12b. The process must be adjusted for 
each pattern; the sequence of photography, developing, and 
the printing of the positive image is highly dependent on the 
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ability of the photographer. Concurrently, the quality of 
reproduction of these patteiTis in printed form depends on th> 
plate preparation, the printing press, and the paper used. 

B. Analysis. 

A 30 X 30 square grid is placed over the centra) 3*0 x 3-0 
unit square of the raw pattern by program DATAGRID. 

The identification number of an atom and the impact point and 
impact area of the run in which that atom was sputtered ai>e 
recorded for each atom in the grid square. The total number 
of atoms, the total energy, and the average energy per atom 
for each grid square are printed in separate square arrays. 
The individual grid square data provides for correlation 
between an atom's crystal location and its deposit point in 
the pattern; the square arrays of the number of atoms and 
the total and average energy densities may be compared 
directly with the point plot or smooth pattern for analysis ol' 
pattern features. 
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FPTF(X )=EXPF( AC+EX8*X ) 
PFPTF( X )=EXPF (PAC+PEXB*X) 
POTF(X ) =EXPF( EXA+EXB*X ) 
PPOTF( X )=EXPF ( PEXA+PEXB*X I 
DO 2 1=1,2000 



9010 FORMAK 10A8 ) 
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THIS IS A LATTICE GENERATOR FOR THE (111) ORIENTATION. 

THE CRYSTAL IS DEVELOPED IN THE ORDER, Z FOLLOWED BY Y, FOLLOWED BY 
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FUNCTION TIMO(IT) 

XlT = n COG 

TIM0=XIT/60. 
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APPENDIX D 



Section 1: Memory block allocations and functions are 

established. All storage cells are zeroed. 

Input and output formats are specified. 

Constants are set, and the target material, 
ion species, and crystal face to be sputtered 
are read in. 

Section 2: The ion energy, impact area, and impact. 

point are read in. Constants peculiar to the 
run are established. The appropriate sub- 
routine is called to build the crystal, and 
crystal boundaries are calculated. Subroutine 
START is called to position the ion. The 
initial value of ^t is calculated. Initial 
coordinates of all atoms are assigned, and 
their velocity components are zeroed (except 
for the ion). The initial coordinates of the 
crystal atoms are printed. 

Section 3 ’ Forces are calculated by calling subroutine 

STEP. Atoms are moved to their intermediate 
positions in the two step cycle. Subroutine 
STEP is called again, and atoms are moved 
to their final positions. Final velocities are 
then computed. Force components are 



zeroed in preparation for the next timestep. 



The maximum kinetic energy is determined for 
the calculation of^t. 

Section l+J The time remaining until cutoff is determined. 

If there is insufficient time to complete 
another timestep, terminal data is printed. 
Potential energy is calculated and summed 
with kinetic energy to give the total energy 
for energy balance check (manual check). The 
data for atoms having potential energy greater 
than the thermal energy are printed. 

Section 5* Atoms which have kinetic energy but not 

potential energy are assumed to be free of 
the crystal. They are assigned LCUT=1. 

The surface through which an atom exited is 
determined and a code assigned. A maximum 
potential energy is found (for Lt calculation) 
among atoms which do not have LCUT=1. If 
the maximum potential energy is less than a 
minimum value the terminating process begins. 

If this energy is greater than the minimum 
value, another timestep begins. 

Section 6: Pertinent data for all atoms is printed. Atoms 

which have exited through the front of the 
crystal or will exit through the front are 
assumed to be sputtered. Data for these 
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atoms arepi"inted and a data card punched 



for each atom. A new data card is then 
read into initiate another run. 
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APPENDIX E 





Glossary for FCCSPUT 


AC 


Parameter for target force function correction. 


AIX 

AIY 

AIS 


Floating point form of IX, lY, IE. 


ALFA 


Cosine of the angle between vectors’ Rl, R2. 


ALFA2 


ALFA squared. 


B 


Reciprocal of magnitude of atom velocity. 


BENGY 


Energy which an atom within the crystal at shut- 
down must have to be considered sputtered. 


BULLET 


Variable representing primary material. 


BX 

BS 


Unsealed x, z coordinates of the impact point 


CELS 


Frictional force multiplier. (See CVS) 


COX 

COY 

COE 


Direction cosines of primary velocity vector. 


COXI 

COYI 

COSI 


Negative values of COX, COY, COE. 


CVB 


A constant. 


CVD 

CVE 

CVM 


Converts meters to angstrom units. 
Converts electron volts to joules. 
Converts atomic mass units to kilograms. 


CVR 


Converts lattice units to angstrom units. 


DCOX 

DCOY 

DC02 


Direction cosines of sputtered atom velocity vector. 


DFF 

DIST 


Distance difference between nearest neighbor 
distance and actual atom ,^diff erences . 
Distance between any two atoms. 
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DRX 

DRY 

DRS 


X, y, z components of DIST. 


DT 


Length of timestep in seconds. 


DTI 


Number of lattice units most energetic atom 
may move in one timestep. 


DTOD 


DT/CVD - a ratio used to avoid repeated division. 


DTOM 


DT/PTMAS - a ratio used to avoid repeated 
division . 


DTOMB 


DT/PGMAS - a ratio used to avoid repeated 
division . 


DX 

DY 

DS 


X, y, z distances atom has moved from initial 
position . 


ECUT 


A lower limit on an atom’s potential energy. If 
energy is less than or equal to ECUT the program 
shuts down. 


El 


A cutoff energy. 


EMAX 


The maximum energy encountered in any cycle - 


EV 


Primary energy in electron volts. 


EVR 


Primary energy in kilo-electron volts. 


EXA 

EXB 


Potential function parameters. 


FA 


The component force increment on an atom. 


FAC 


The minimum distance the primary is positioned 
in front of the first xz plane at start time. 


FM 


A small number used in checking potential energy 
zero point. 


FOD 


FORCE/DIST - a ratio used to avoid repeated 
division. 


FORCE 


Numerical value of the force function with a 
variable parameter. 
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FORF 


Target atom force function. 


FPTC 


The corrective force value at ROE. 


FPTF 


The corrective force function. 


FRC 


Numerical value of the target force function at 
ROE. 


FX 

FY 

FS 


X, y, z components of total force on an atom. 


FXA- 


Force function parameter. 


GMAS 


Target atom mass (in a.m.u.) 


HDTOD 


1/2 DTOD - a ratio used to avoid repeated 
division . 


HDTOM 


1/2 DTOM — a ratio used to avoid repeated 
division. 


HDTOMB 


1/2 DTOMB — a ratio used to avoid repeated 
division. 


HGMAS 


1/2 GMAS - a ratio used to avoid repeated 
division . 


HTMAS 


1/2 TMAS - a ratio used to avoid repeated 
division. 


ICUT 


Used to provide output prior to time limit 
shutdown. 


IDUM 


Dummy variable. 


IHB 
IHS . 
IHT 
IHl 


Alpha-numeric arrays for titling. 


IN 


Odd-even integer used to determine atom site 
establishment . 


INDEX 


Integer (0 or 1 ) used in determining dynamics 
cycle step. 


I NOW 


Time program has been running in seconds. 
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IP 


Subscript value of atom. Used in subroutine 

STEP. 


ISHUT 


Time left prior to time limit. 


ISPX 

ISP2 


Fixed poinr values of SPX, SPS. 


ISS 


Subscript value of most energetic atom. 


IT 


Unsealed fixed point x coordinate used in lattice 
generation. Also a dummy variable in function 
TIMEF. 


ITT 


Odd-even integer used to determine atom site 
establishment . 


IX 

lY 

IZ 


Number of x, y, z planes of crystal. 


IXP 

lYP 

ISP 


Crystal dimensions in x, y, z. 


aSHUT 


Cutoff variable based on total potential energy 
of crystal. 


JT 


Unsealed y coordinate used in crystal generation. 


aTS 

JTT 


Variables used to establish atom sites. 


KCUT 


Identifies exit point of atom. 


KT 


Unsealed z coordinate used to establish atom site 


LCUT 


Used to identify atoms which are not included in 
calculations . 


LL 


The highest numbered atom in the crystal. 


LS 


Sum of the Miller index integers. 


LSS 


Used to identify type of surface, i.e., regular, 
stub, vacancy. 


M 


An integer used to begin atom numbering. 
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N 


Subscript of the atom to be removed for 
vacancy surface. 


ND 


Data output increment. 


NPAGE 


Page numbering variable. 


NS 


Initial print statement cycle . 


NSHUT 


Cutoff variable based on too long a timestep 


NT 


Timestep . 


NTT 


Timestep limit before shutdown. 


PAG 


Same as AC except applicable to primary. 


PEXA 

PEXB 


Primary force function parameters. 


PFIV 


A constant = 0.5* 


PFORF 


Primary force function. 


PFPTC 


PFPTF evaluated at ROE. 


PFPTF 


Primary corrective force function. 


PFRC 


PFORF evaluated at ROE. 


PFXA 


Primary force function parameter. 


PGMAS 


Primary mass in kilograms. 


PKE 


Kinetic energy of an atom. 


PKEY 


Y component of kinetic energy of an atom. 


pea’ 


Crystal plane (alphanumeric variable). 


PLANE 


Same as PLA. 


PNUM 


Impact point (alphanumeric variable). 


POT 


Potential energy between two atoms. 


POTF 


Target potential function. 


PPE 


Potential energy of an atom. 
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PPOTF 


Primary potential function. 


PPTC 


PPOTF evaluated at ROE. 


PRI 


Chemical symbol for primary material. 


PTC 


POTF evaluated at ROE. 


PTE 


Total energy of an atom (potential t kinetic). 


PTMAS 


Target mass in kilograms. 


QM 


A small number used in checking kinetic energy 
zero point. 


QUIT 


Cutoff variable checked against total potential 
energy . 


RBX 

RBZ 


Unsealed x, z coordinates of impact area 
reference point 


ROE 


Nearest neighbor distance. 


ROEM 


ROE ~ DTI (one timestep distance less than 
n . n . distance ) . 


ROE2 


ROE squared. 


RX 

RY 

RE 


X, y, z coordinates of atom at any time. 


RXBND 

RYBND 

REBND 


X, y, z coordinates of crystal boundaries other 
than zero. 


RXI 

RYI 

R2I 


X, y, z coordinates of atom's initial position. 


RXK 

RYK 

REK 


X, y, z coordinates of temporary position of 
atom during force cycle. 


RXS 

RYS 

RES 


X, y, z coordinates of impact point. 


R1 


Vector from impact point to initial primary position 
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RIX 

RIY 

R13 


X, y, z coordinates of initial primary position. 


R1R2 


Scalar product of vectors Rl, R2. 


R2 


Magnitude of vector from impact point to first 
atom hit by primary. 


R2SQ 


R2 squared. 


R3 


Magnitude of vector from impact point of 
primary start position. 


SAVE 


1/2 POT. 


sex 

SCY 

SC3 


X, y, z coordinate scale factors. 


SLOW 


Cutoff variable checked against a long DT . 


SPX 

SPS 


X, z distance from impact area reference point 
to impact point. 


SSCS 


A z scale factor used for (ill) plane lattice 
generation . 


SUR 


Plane (alphanumeric variable). 


TAR 


Chemical symbol for target material. 


TARGET 


Target material (alphanumeric variable)* 


TE 


Total energy of crystal atoms (kinetic + 
potential) . 


TEMP 


Temperature of lattice in degrees Kelvin. 


THERM 


Thermal energy of atom. 


TI 


Computer time program has been running. 


TIME 


Elapsed problem time. 


TIMO 


A function to convert seconds to minutes. 


TMAS 


Target atom mass in kilograms. 
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TPKE 


Total kinetic energy of crystal atoms. 


TPOT 


Total potential energy of crystal atoms. 


VOL 


Magnitude of primary velocity vector. 


VSS 


Storage variable for velocity components. 


VX 

VY 

VS 


X, y, z components of atom velocity. 


X 


Unsealed x coordinate used in crystal generation. 


XIT 


Floating point form of IT used in function TIMO 


Y 


Unsealed y coordinate used in crystal generation. 


Z 


Unsealed z coordinate used in crystal generation. 


SE 


A constant =-• 0.0 


SP 


Floating point form of <JTT. 
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BO 110 I=2»3000 

READ 5*PNUM»NR,VX»VY»VZ tPKEY 

IF(NR) 100,120,100 

100 WRITE TAPE 2 , PNUM , NR , VX , VY , VZ , PKE Y 
no CONTINUE 
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PKE=(VX*VX + VY*VY + VZ*VZ)*HTMAS 
II=PKE + 0.49 
IF ( 1 1-500 ) 3060*305 0.3050 
3050 11=500 

PRINT 10,PN0M»NR»Pf;EY 
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PROGRAM DATAGRID 
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PRINT ^4l*PLANE 
PRINT 42*EVK,EGYCUT 
IF(C0y,EQ.1.0)211,212 
211 PRINT 9U 
GO TO 213 
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Erosion of the force function 
at r - 2r^ results in a loss 
of energy in the model . Thi s 
is corrected in the STEP subroutine. 

Figure 1 
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Figure 44a 
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Figure 44b 
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Figure 46 
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Figure 51c 
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Figure 61b 
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Figure 62a 
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Figure 62b 
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Figure 63a 
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Figure 63b 
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